GGW-mod2: Moral Followup (to Study 1)

Setup

## Warning: package 'ggplot2' was built under R version 3.2.4

conditionsExcluded	outlierHandling	EFA_correlation
none	keep	cor

back to TOC

Data preparation

## Warning in read.table(file = file, header = header, sep = sep, quote
## = quote, : incomplete final line found by readTableHeader on '/Users/
## kweisman/Documents/Research (Stanford)/Projects/GGW-mod/ggw-mod2_moral/
## mturk/v1 follow up 1 (robot)/Batch_2292738_batch_results.csv'

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning in ifelse(as.numeric(as.character(structure(c(48L, 47L, 50L, 43L, :
## NAs introduced by coercion

## Warning in ifelse(as.numeric(as.character(structure(c(48L, 47L, 50L, 43L, :
## NAs introduced by coercion

## Warning in ifelse(as.numeric(as.character(structure(c(48L, 47L, 50L, 43L, :
## NAs introduced by coercion

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining factors with different levels, coercing to character
## vector

## Warning: joining character vector and factor, coercing into character
## vector

## Warning: joining factor and character vector, coercing into character
## vector

back to TOC

Study 1 (2015-12-15, 2 conditions, between-subjects)

Demographics

Study 1: Sample size

condition	n
beetle	200
robot	205
all	405

Study 1: Study duration (minutes)

condition	min_duration	max_duration	median_duration	mean_duration	sd_duration
beetle	0.77	119.12	2.77	3.92	8.50
robot	0.83	15.77	2.58	3.24	2.11
all	0.77	119.12	2.67	3.58	6.16

## Joining by: c("condition", "min_age", "max_age", "median_age", "mean_age", "sd_age")

## Warning: joining character vector and factor, coercing into character
## vector

Study 1: Participant age (years; approximate)

condition	min_age	max_age	median_age	mean_age	sd_age
beetle	19	70	31	34.24	10.64
robot	20	75	32	35.11	11.20
all	19	75	32	34.68	10.92

Study 1: Participant gender

condition	m	f	other
beetle	119	81	0
robot	121	82	2
all	240	163	2

Study 1: Participant race/ethnicity

condition	asian_east	asian_other	asian_south	black	hispanic	middle_eastern	multiracial	native_american	other_prefno	pac_islander	white
beetle	20	0	6	8	9	1	11	1	1	1	142
robot	16	2	2	15	9	1	9	1	1	1	148
all	36	2	8	23	18	2	20	2	2	2	290

Study 1: Participant religion

condition	buddhism	christianity	hinduism	islam	judaism	multireligious	none	other	other_prefno
beetle	3	81	2	2	4	0	101	3	4
robot	2	87	0	0	5	3	100	3	5
all	5	168	2	2	9	3	201	6	9

back to TOC

Exploratory factor analysis

Run EFA with varimax rotation

## Factor Analysis using method =  minres
## Call: fa(r = d1_all, nfactors = 3, rotate = "varimax", scores = "regression", 
##     fm = "minres", cor = chosenCorType)
## Standardized loadings (pattern matrix) based upon correlation matrix
##                  MR1   MR2   MR3   h2   u2 com
## happy           0.55  0.65  0.10 0.74 0.26 2.0
## depressed       0.37  0.78  0.03 0.75 0.25 1.4
## fear            0.82  0.28  0.06 0.76 0.24 1.2
## angry           0.58  0.57  0.08 0.67 0.33 2.0
## calm            0.65  0.41  0.17 0.62 0.38 1.9
## sounds          0.06 -0.05  0.61 0.38 0.62 1.0
## seeing          0.36 -0.07  0.61 0.51 0.49 1.7
## temperature     0.21 -0.06  0.66 0.48 0.52 1.2
## odors           0.45 -0.01  0.43 0.39 0.61 2.0
## depth           0.11  0.11  0.62 0.40 0.60 1.1
## computations   -0.74  0.19  0.44 0.78 0.22 1.8
## thoughts        0.50  0.55  0.22 0.61 0.39 2.3
## reasoning      -0.06  0.47  0.57 0.55 0.45 2.0
## remembering    -0.20  0.17  0.72 0.59 0.41 1.3
## beliefs         0.11  0.76  0.12 0.60 0.40 1.1
## hungry          0.93  0.01 -0.08 0.87 0.13 1.0
## tired           0.83  0.23  0.10 0.76 0.24 1.2
## pain            0.93  0.10  0.01 0.87 0.13 1.0
## nauseated       0.65  0.50  0.09 0.68 0.32 1.9
## safe            0.70  0.36  0.13 0.63 0.37 1.6
## love            0.37  0.81  0.06 0.79 0.21 1.4
## recognizing    -0.29  0.29  0.71 0.67 0.33 1.7
## communicating  -0.02  0.20  0.65 0.46 0.54 1.2
## guilt           0.26  0.80  0.02 0.70 0.30 1.2
## disrespected    0.25  0.78  0.04 0.68 0.32 1.2
## free_will       0.70  0.37  0.09 0.63 0.37 1.5
## choices         0.25  0.18  0.60 0.46 0.54 1.5
## self_restraint  0.24  0.55  0.31 0.45 0.55 2.0
## intentions      0.54  0.35  0.27 0.48 0.52 2.3
## goal            0.09  0.17  0.62 0.42 0.58 1.2
## conscious       0.70  0.36  0.12 0.64 0.36 1.6
## self_aware      0.52  0.48  0.22 0.55 0.45 2.3
## desires         0.69  0.40  0.11 0.66 0.34 1.7
## embarrassed     0.19  0.85 -0.01 0.76 0.24 1.1
## emo_recog       0.06  0.70  0.29 0.58 0.42 1.3
## joy             0.51  0.70  0.10 0.76 0.24 1.9
## morality       -0.04  0.60  0.32 0.47 0.53 1.5
## personality     0.23  0.66  0.31 0.58 0.42 1.7
## pleasure        0.74  0.43  0.11 0.74 0.26 1.7
## pride           0.28  0.85  0.05 0.80 0.20 1.2
## 
##                        MR1  MR2  MR3
## SS loadings           9.92 9.81 5.17
## Proportion Var        0.25 0.25 0.13
## Cumulative Var        0.25 0.49 0.62
## Proportion Explained  0.40 0.39 0.21
## Cumulative Proportion 0.40 0.79 1.00
## 
## Mean item complexity =  1.6
## Test of the hypothesis that 3 factors are sufficient.
## 
## The degrees of freedom for the null model are  780  and the objective function was  36.39 with Chi Square of  14185
## The degrees of freedom for the model are 663  and the objective function was  3.81 
## 
## The root mean square of the residuals (RMSR) is  0.03 
## The df corrected root mean square of the residuals is  0.03 
## 
## The harmonic number of observations is  405 with the empirical chi square  500.89  with prob <  1 
## The total number of observations was  405  with MLE Chi Square =  1475.92  with prob <  8.7e-64 
## 
## Tucker Lewis Index of factoring reliability =  0.928
## RMSEA index =  0.057  and the 90 % confidence intervals are  0.051 0.059
## BIC =  -2504.66
## Fit based upon off diagonal values = 1
## Measures of factor score adequacy             
##                                                 MR1  MR2  MR3
## Correlation of scores with factors             0.98 0.98 0.95
## Multiple R square of scores with factors       0.97 0.95 0.91
## Minimum correlation of possible factor scores  0.93 0.91 0.81

## Joining by: "subid"

## Warning: joining factor and character vector, coercing into character
## vector

## Joining by: "mturkcode"
## Joining by: "workerID"

## Warning: joining factor and character vector, coercing into character
## vector

back to TOC

Study 1 Moral Follow-up (robot: 2016-02-26; beetle: 2016-03-01)

Data preparation

Demographics

Study 1: Sample size

condition	n
beetle	111
robot	121
all	232

Study 1: Study duration (minutes)

condition	min_duration	max_duration	median_duration	mean_duration	sd_duration
beetle	0.77	23.67	2.68	3.36	3.22
robot	0.82	11.45	2.48	3.10	1.95

back to TOC

Correlations with general moral concern

## 
##  Pearson's product-moment correlation
## 
## data:  MR1 and concern_general
## t = 4.4105, df = 229, p-value = 1.587e-05
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  0.156382 0.394636
## sample estimates:
##       cor 
## 0.2798115

## 
##  Pearson's product-moment correlation
## 
## data:  MR2 and concern_general
## t = 3.4838, df = 229, p-value = 0.0005921
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  0.09811329 0.34348073
## sample estimates:
##       cor 
## 0.2243498

## 
##  Pearson's product-moment correlation
## 
## data:  MR3 and concern_general
## t = 0.047144, df = 229, p-value = 0.9624
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.1260130  0.1321399
## sample estimates:
##         cor 
## 0.003115362

## 
##  Pearson's product-moment correlation
## 
## data:  MR1 and MR2
## t = -0.13729, df = 229, p-value = 0.8909
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.1379882  0.1201462
## sample estimates:
##          cor 
## -0.009072134

## 
##  Pearson's product-moment correlation
## 
## data:  MR1 and MR3
## t = 0.48221, df = 229, p-value = 0.6301
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.09762977  0.16026807
## sample estimates:
##        cor 
## 0.03184926

## 
##  Pearson's product-moment correlation
## 
## data:  MR2 and MR3
## t = 0.040367, df = 229, p-value = 0.9678
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.1264537  0.1316998
## sample estimates:
##       cor 
## 0.0026675

## Call: paired.r(xy = rp_MR1.gen$estimate, xz = rp_MR2.gen$estimate, 
##     yz = rp_MR1.MR2$estimate, n = d_moral_merged_n)
## [1] "test of difference between two correlated  correlations"
## t = 0.62  With probability =  0.53

## Call: paired.r(xy = rp_MR1.gen$estimate, xz = rp_MR3.gen$estimate, 
##     yz = rp_MR1.MR3$estimate, n = d_moral_merged_n)
## [1] "test of difference between two correlated  correlations"
## t = 3.12  With probability =  0

## Call: paired.r(xy = rp_MR2.gen$estimate, xz = rp_MR3.gen$estimate, 
##     yz = rp_MR2.MR3$estimate, n = d_moral_merged_n)
## [1] "test of difference between two correlated  correlations"
## t = 2.42  With probability =  0.02

## Warning in cor.test.default(MR1, concern_general, method = "spearman"):
## Cannot compute exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  MR1 and concern_general
## S = 1435500, p-value = 3.123e-06
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##       rho 
## 0.3012511

## Warning in cor.test.default(MR2, concern_general, method = "spearman"):
## Cannot compute exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  MR2 and concern_general
## S = 1682900, p-value = 0.005847
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##       rho 
## 0.1808292

## Warning in cor.test.default(MR3, concern_general, method = "spearman"):
## Cannot compute exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  MR3 and concern_general
## S = 2162500, p-value = 0.4261
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##         rho 
## -0.05261743

## Warning in cor.test.default(MR1, MR2, method = "spearman"): Cannot compute
## exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  MR1 and MR2
## S = 2476500, p-value = 0.001691
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##        rho 
## -0.2054888

## Warning in cor.test.default(MR1, MR3, method = "spearman"): Cannot compute
## exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  MR1 and MR3
## S = 2115200, p-value = 0.6541
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##        rho 
## -0.0296376

## Warning in cor.test.default(MR2, MR3, method = "spearman"): Cannot compute
## exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  MR2 and MR3
## S = 2106800, p-value = 0.6994
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##        rho 
## -0.0255419

## Call: paired.r(xy = rs_MR1.genrank$estimate, xz = rs_MR2.genrank$estimate, 
##     yz = rs_MR1.MR2$estimate, n = d_moral_merged_n)
## [1] "test of difference between two correlated  correlations"
## t = 1.24  With probability =  0.22

## Call: paired.r(xy = rs_MR1.genrank$estimate, xz = rs_MR3.genrank$estimate, 
##     yz = rs_MR1.MR3$estimate, n = d_moral_merged_n)
## [1] "test of difference between two correlated  correlations"
## t = 3.9  With probability =  0

## Call: paired.r(xy = rs_MR2.genrank$estimate, xz = rs_MR3.genrank$estimate, 
##     yz = rs_MR2.MR3$estimate, n = d_moral_merged_n)
## [1] "test of difference between two correlated  correlations"
## t = 2.51  With probability =  0.01

## 
##  Pearson's product-moment correlation
## 
## data:  MR1 and concern_general
## t = 3.1616, df = 119, p-value = 0.001991
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  0.1050907 0.4352364
## sample estimates:
##       cor 
## 0.2783656

## 
##  Pearson's product-moment correlation
## 
## data:  MR2 and concern_general
## t = 3.042, df = 119, p-value = 0.002892
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  0.09464914 0.42664831
## sample estimates:
##       cor 
## 0.2686077

## 
##  Pearson's product-moment correlation
## 
## data:  MR3 and concern_general
## t = 0.0061249, df = 119, p-value = 0.9951
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.1779529  0.1790401
## sample estimates:
##          cor 
## 0.0005614675

## 
##  Pearson's product-moment correlation
## 
## data:  MR1 and MR2
## t = 13.648, df = 119, p-value < 2.2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  0.7002746 0.8422012
## sample estimates:
##       cor 
## 0.7811323

## 
##  Pearson's product-moment correlation
## 
## data:  MR1 and MR3
## t = 0.22997, df = 119, p-value = 0.8185
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.1580147  0.1988248
## sample estimates:
##       cor 
## 0.0210763

## 
##  Pearson's product-moment correlation
## 
## data:  MR2 and MR3
## t = -0.52753, df = 119, p-value = 0.5988
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.2248602  0.1313264
## sample estimates:
##         cor 
## -0.04830237

## Call: paired.r(xy = rp_ROBOT_MR1.gen$estimate, xz = rp_ROBOT_MR2.gen$estimate, 
##     yz = rp_ROBOT_MR1.MR2$estimate, n = d_moral_merged_n)
## [1] "test of difference between two correlated  correlations"
## t = 0.23  With probability =  0.82

## Call: paired.r(xy = rp_ROBOT_MR1.gen$estimate, xz = rp_ROBOT_MR3.gen$estimate, 
##     yz = rp_ROBOT_MR1.MR3$estimate, n = d_moral_merged_n)
## [1] "test of difference between two correlated  correlations"
## t = 3.11  With probability =  0

## Call: paired.r(xy = rp_ROBOT_MR2.gen$estimate, xz = rp_ROBOT_MR3.gen$estimate, 
##     yz = rp_ROBOT_MR2.MR3$estimate, n = d_moral_merged_n)
## [1] "test of difference between two correlated  correlations"
## t = 2.89  With probability =  0

## Warning in cor.test.default(MR1, concern_general, method = "spearman"):
## Cannot compute exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  MR1 and concern_general
## S = 229350, p-value = 0.01387
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##       rho 
## 0.2231781

## Warning in cor.test.default(MR2, concern_general, method = "spearman"):
## Cannot compute exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  MR2 and concern_general
## S = 209340, p-value = 0.001205
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##       rho 
## 0.2909529

## Warning in cor.test.default(MR3, concern_general, method = "spearman"):
## Cannot compute exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  MR3 and concern_general
## S = 311950, p-value = 0.5376
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##         rho 
## -0.05658593

## Warning in cor.test.default(MR1, MR2, method = "spearman"): Cannot compute
## exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  MR1 and MR2
## S = 231830, p-value = 0.018
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##       rho 
## 0.2147718

## Warning in cor.test.default(MR1, MR3, method = "spearman"): Cannot compute
## exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  MR1 and MR3
## S = 283930, p-value = 0.6765
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##        rho 
## 0.03831133

## Warning in cor.test.default(MR2, MR3, method = "spearman"): Cannot compute
## exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  MR2 and MR3
## S = 360410, p-value = 0.01498
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##       rho 
## -0.220733

## Call: paired.r(xy = rs_ROBOT_MR1.genrank$estimate, xz = rs_ROBOT_MR2.genrank$estimate, 
##     yz = rs_ROBOT_MR1.MR2$estimate, n = d_moral_merged_n)
## [1] "test of difference between two correlated  correlations"
## t = -0.86  With probability =  0.39

## Call: paired.r(xy = rs_ROBOT_MR1.genrank$estimate, xz = rs_ROBOT_MR3.genrank$estimate, 
##     yz = rs_ROBOT_MR1.MR3$estimate, n = d_moral_merged_n)
## [1] "test of difference between two correlated  correlations"
## t = 3.13  With probability =  0

## Call: paired.r(xy = rs_ROBOT_MR2.genrank$estimate, xz = rs_ROBOT_MR3.genrank$estimate, 
##     yz = rs_ROBOT_MR2.MR3$estimate, n = d_moral_merged_n)
## [1] "test of difference between two correlated  correlations"
## t = 3.49  With probability =  0

## 
##  Pearson's product-moment correlation
## 
## data:  MR1 and concern_general
## t = 1.1775, df = 108, p-value = 0.2416
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.07626003  0.29364023
## sample estimates:
##       cor 
## 0.1125891

## 
##  Pearson's product-moment correlation
## 
## data:  MR2 and concern_general
## t = 2.4254, df = 108, p-value = 0.01695
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  0.04181352 0.39759715
## sample estimates:
##       cor 
## 0.2272756

## 
##  Pearson's product-moment correlation
## 
## data:  MR3 and concern_general
## t = 1.5122, df = 108, p-value = 0.1334
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.04444539  0.32253964
## sample estimates:
##       cor 
## 0.1439944

## 
##  Pearson's product-moment correlation
## 
## data:  MR1 and MR2
## t = -1.7993, df = 108, p-value = 0.07477
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.34676327  0.01719213
## sample estimates:
##        cor 
## -0.1705985

## 
##  Pearson's product-moment correlation
## 
## data:  MR1 and MR3
## t = 17.792, df = 108, p-value < 2.2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  0.8066654 0.9044886
## sample estimates:
##      cor 
## 0.863485

## 
##  Pearson's product-moment correlation
## 
## data:  MR2 and MR3
## t = 0.28544, df = 108, p-value = 0.7759
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.1606108  0.2135997
## sample estimates:
##        cor 
## 0.02745637

## Call: paired.r(xy = rp_BEETLE_MR1.gen$estimate, xz = rp_BEETLE_MR2.gen$estimate, 
##     yz = rp_BEETLE_MR1.MR2$estimate, n = d_moral_merged_n)
## [1] "test of difference between two correlated  correlations"
## t = -1.16  With probability =  0.25

## Call: paired.r(xy = rp_BEETLE_MR1.gen$estimate, xz = rp_BEETLE_MR3.gen$estimate, 
##     yz = rp_BEETLE_MR1.MR3$estimate, n = d_moral_merged_n)
## [1] "test of difference between two correlated  correlations"
## t = -0.92  With probability =  0.36

## Call: paired.r(xy = rp_BEETLE_MR2.gen$estimate, xz = rp_BEETLE_MR3.gen$estimate, 
##     yz = rp_BEETLE_MR2.MR3$estimate, n = d_moral_merged_n)
## [1] "test of difference between two correlated  correlations"
## t = 0.93  With probability =  0.35

## Warning in cor.test.default(MR1, concern_general, method = "spearman"):
## Cannot compute exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  MR1 and concern_general
## S = 199790, p-value = 0.302
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##        rho 
## 0.09930608

## Warning in cor.test.default(MR2, concern_general, method = "spearman"):
## Cannot compute exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  MR2 and concern_general
## S = 177680, p-value = 0.03715
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##       rho 
## 0.1989936

## Warning in cor.test.default(MR3, concern_general, method = "spearman"):
## Cannot compute exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  MR3 and concern_general
## S = 2e+05, p-value = 0.3066
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##        rho 
## 0.09836269

## Warning in cor.test.default(MR1, MR2, method = "spearman"): Cannot compute
## exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  MR1 and MR2
## S = 286800, p-value = 0.001898
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##        rho 
## -0.2929481

## Warning in cor.test.default(MR1, MR3, method = "spearman"): Cannot compute
## exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  MR1 and MR3
## S = 35168, p-value < 2.2e-16
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##       rho 
## 0.8414527

## Warning in cor.test.default(MR2, MR3, method = "spearman"): Cannot compute
## exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  MR2 and MR3
## S = 240360, p-value = 0.3852
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##         rho 
## -0.08360158

## Call: paired.r(xy = rs_BEETLE_MR1.genrank$estimate, xz = rs_BEETLE_MR2.genrank$estimate, 
##     yz = rs_BEETLE_MR1.MR2$estimate, n = d_moral_merged_n)
## [1] "test of difference between two correlated  correlations"
## t = -0.96  With probability =  0.34

## Call: paired.r(xy = rs_BEETLE_MR1.genrank$estimate, xz = rs_BEETLE_MR3.genrank$estimate, 
##     yz = rs_BEETLE_MR1.MR3$estimate, n = d_moral_merged_n)
## [1] "test of difference between two correlated  correlations"
## t = 0.03  With probability =  0.98

## Call: paired.r(xy = rs_BEETLE_MR2.genrank$estimate, xz = rs_BEETLE_MR3.genrank$estimate, 
##     yz = rs_BEETLE_MR2.MR3$estimate, n = d_moral_merged_n)
## [1] "test of difference between two correlated  correlations"
## t = 1.06  With probability =  0.29

Exploratory plots

Plots with general moral concern

## Warning: Removed 1 rows containing non-finite values (stat_smooth).

## Warning: Removed 1 rows containing missing values (geom_point).

## Warning: Removed 1 rows containing non-finite values (stat_smooth).

## Warning: Removed 1 rows containing missing values (geom_point).

## Warning: Removed 1 rows containing non-finite values (stat_smooth).

## Warning: Removed 1 rows containing missing values (geom_point).

## Warning: Removed 3 rows containing non-finite values (stat_smooth).

## Warning: Removed 3 rows containing missing values (geom_point).

## Warning: Removed 1 rows containing non-finite values (stat_bin).

## Warning: Removed 1 rows containing missing values (geom_bar).

## Warning: Removed 1 rows containing missing values (geom_bar).

## Warning: Removed 1 rows containing non-finite values (stat_bin).

## Warning: Removed 1 rows containing missing values (geom_bar).

## Warning: Removed 1 rows containing missing values (geom_bar).

## Warning: Removed 1 rows containing non-finite values (stat_bin).

## Warning: Removed 1 rows containing missing values (geom_bar).

## Warning: Removed 1 rows containing missing values (geom_bar).

Regressions with general moral concern

## Analysis of Variance Table
## 
## Model 1: concern_general ~ MR1 + MR2 + MR3 + condition
## Model 2: concern_general ~ (MR1 + MR2 + MR3) * condition
## Model 3: concern_general ~ (MR1 + MR2 + MR3 + condition)^2
## Model 4: concern_general ~ (MR1 + MR2 + MR3)^2 * condition
## Model 5: concern_general ~ (MR1 + MR2 + MR3)^3 * condition
## Model 6: concern_general ~ (MR1 + MR2 + MR3 + condition)^3
##   Res.Df    RSS Df Sum of Sq      F   Pr(>F)   
## 1    226 111315                                
## 2    223 111154  3     161.2 0.1135 0.952138   
## 3    220 104031  3    7123.2 5.0166 0.002213 **
## 4    217 103464  3     567.1 0.3994 0.753576   
## 5    215 101761  2    1702.3 1.7983 0.168065   
## 6    216 102059 -1    -297.6 0.6288 0.428680   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## lm(formula = concern_general ~ (MR1 + MR2 + MR3 + condition)^2, 
##     data = d_moral_merged)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -49.836 -14.847  -4.775  11.270  68.561 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         34.9356     6.5930   5.299 2.82e-07 ***
## MR1                 17.7455     6.7744   2.620 0.009419 ** 
## MR2                  5.9861     2.5214   2.374 0.018448 *  
## MR3                 -0.1649     2.7975  -0.059 0.953046    
## conditionrobot      11.3630     6.5462   1.736 0.083995 .  
## MR1:MR2            -12.1356     3.3463  -3.627 0.000357 ***
## MR1:MR3              1.0616     2.9577   0.359 0.719994    
## MR1:conditionrobot  12.4343     6.9827   1.781 0.076336 .  
## MR2:MR3              1.8989     2.7016   0.703 0.482873    
## MR2:conditionrobot -10.2120     3.7136  -2.750 0.006457 ** 
## MR3:conditionrobot   1.3525     3.8584   0.351 0.726268    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 21.75 on 220 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.1868, Adjusted R-squared:  0.1499 
## F-statistic: 5.055 on 10 and 220 DF,  p-value: 1.26e-06

## Analysis of Variance Table
## 
## Model 1: concern_general ~ (MR1 + MR2 + condition)^2
## Model 2: concern_general ~ (MR1 + MR2)^2 * condition
## Model 3: concern_general ~ (MR1 + MR2 + condition)^3
##   Res.Df    RSS Df Sum of Sq      F Pr(>F)
## 1    224 104328                           
## 2    223 103786  1    541.91 1.1644 0.2817
## 3    223 103786  0      0.00

## Analysis of Variance Table
## 
## Model 1: concern_general ~ (MR1 + MR2 + condition)^2
## Model 2: concern_general ~ (MR1 + MR2 + MR3 + condition)^2
##   Res.Df    RSS Df Sum of Sq      F Pr(>F)
## 1    224 104328                           
## 2    220 104031  4    296.69 0.1569 0.9597

## 
## Call:
## lm(formula = concern_general ~ (MR1 + MR2 + condition)^2, data = d_moral_merged)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -47.813 -14.897  -4.131  11.159  68.378 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)          35.224      5.247   6.713 1.55e-10 ***
## MR1                  16.147      5.840   2.765 0.006167 ** 
## MR2                   6.077      2.450   2.480 0.013862 *  
## conditionrobot        9.455      5.241   1.804 0.072565 .  
## MR1:MR2             -11.004      2.874  -3.829 0.000167 ***
## MR1:conditionrobot   12.299      6.011   2.046 0.041907 *  
## MR2:conditionrobot   -8.899      3.203  -2.779 0.005921 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 21.58 on 224 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.1845, Adjusted R-squared:  0.1627 
## F-statistic: 8.447 on 6 and 224 DF,  p-value: 2.84e-08

## Warning: Removed 1 rows containing non-finite values (stat_smooth).

## Warning: Removed 1 rows containing missing values (geom_point).

## Warning: Removed 1 rows containing non-finite values (stat_smooth).

## Warning: Removed 1 rows containing missing values (geom_point).

## Analysis of Variance Table
## 
## Model 1: concern_general ~ rank(MR1) + rank(MR2) + rank(MR3) + condition
## Model 2: concern_general ~ (rank(MR1) + rank(MR2) + rank(MR3)) * condition
## Model 3: concern_general ~ (rank(MR1) + rank(MR2) + rank(MR3) + condition)^2
## Model 4: concern_general ~ (rank(MR1) + rank(MR2) + rank(MR3))^2 * condition
## Model 5: concern_general ~ ((rank(MR1) + rank(MR2) + rank(MR3))^3) * condition
## Model 6: concern_general ~ ((rank(MR1) + rank(MR2) + rank(MR3) + condition)^3)
##   Res.Df    RSS Df Sum of Sq      F  Pr(>F)  
## 1    227 111333                              
## 2    224 110964  3    369.09 0.2534 0.85885  
## 3    221 108699  3   2265.00 1.5552 0.20128  
## 4    218 106584  3   2114.86 1.4521 0.22865  
## 5    216 104864  2   1720.14 1.7716 0.17253  
## 6    217 106369 -1  -1504.51 3.0990 0.07976 .
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## lm(formula = concern_general ~ rank(MR1) + rank(MR2) + rank(MR3) + 
##     condition, data = d_moral_merged)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -42.108 -16.259  -5.906  11.442  74.814 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)     2.953613   5.558039   0.531 0.595652    
## rank(MR1)       0.097392   0.041082   2.371 0.018593 *  
## rank(MR2)       0.082199   0.022388   3.671 0.000301 ***
## rank(MR3)       0.004384   0.024329   0.180 0.857143    
## conditionrobot -1.471971   2.890813  -0.509 0.611114    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 22.15 on 227 degrees of freedom
## Multiple R-squared:  0.1323, Adjusted R-squared:  0.117 
## F-statistic: 8.653 on 4 and 227 DF,  p-value: 1.619e-06

## Analysis of Variance Table
## 
## Model 1: concern_general ~ log(MR1 + 10) + log(MR2 + 10) + log(MR3 + 10) + 
##     condition
## Model 2: concern_general ~ (log(MR1 + 10) + log(MR2 + 10) + log(MR3 + 
##     10)) * condition
## Model 3: concern_general ~ (log(MR1 + 10) + log(MR2 + 10) + log(MR3 + 
##     10) + condition)^2
## Model 4: concern_general ~ (log(MR1 + 10) + log(MR2 + 10) + log(MR3 + 
##     10))^2 * condition
## Model 5: concern_general ~ ((log(MR1 + 10) + log(MR2 + 10) + log(MR3 + 
##     10))^3) * condition
## Model 6: concern_general ~ ((log(MR1 + 10) + log(MR2 + 10) + log(MR3 + 
##     10) + condition)^3)
##   Res.Df    RSS Df Sum of Sq      F   Pr(>F)   
## 1    226 110999                                
## 2    223 110898  3     100.9 0.0709 0.975439   
## 3    220 104136  3    6761.9 4.7508 0.003143 **
## 4    217 103705  3     431.1 0.3029 0.823277   
## 5    215 102004  2    1701.3 1.7930 0.168947   
## 6    216 102349 -1    -345.9 0.7290 0.394157   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## lm(formula = concern_general ~ log(MR1 + 10) + log(MR2 + 10) + 
##     log(MR3 + 10) + condition, data = d_moral_merged)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -42.395 -16.078  -4.117  10.382  73.687 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    -268.4229    75.2895  -3.565 0.000444 ***
## log(MR1 + 10)    67.4818    33.4439   2.018 0.044799 *  
## log(MR2 + 10)    61.0595    16.5333   3.693 0.000278 ***
## log(MR3 + 10)    -0.9552    17.2603  -0.055 0.955918    
## conditionrobot   -0.1739     3.4054  -0.051 0.959320    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 22.16 on 226 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.1324, Adjusted R-squared:  0.117 
## F-statistic:  8.62 on 4 and 226 DF,  p-value: 1.716e-06

## 
## Call:
## lm(formula = concern_general ~ (log(MR1 + 10) + log(MR2 + 10) + 
##     log(MR3 + 10)) * condition, data = d_moral_merged)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -43.800 -16.129  -4.158  10.312  73.740 
## 
## Coefficients:
##                              Estimate Std. Error t value Pr(>|t|)   
## (Intercept)                  -291.035     98.368  -2.959  0.00342 **
## log(MR1 + 10)                  83.143     63.375   1.312  0.19089   
## log(MR2 + 10)                  54.180     27.012   2.006  0.04609 * 
## log(MR3 + 10)                   1.002     28.087   0.036  0.97158   
## conditionrobot                -42.589     98.368  -0.433  0.66547   
## log(MR1 + 10):conditionrobot   24.415     63.375   0.385  0.70042   
## log(MR2 + 10):conditionrobot   -4.795     27.012  -0.178  0.85926   
## log(MR3 + 10):conditionrobot   -0.614     28.087  -0.022  0.98258   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 22.3 on 223 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.1332, Adjusted R-squared:  0.106 
## F-statistic: 4.894 on 7 and 223 DF,  p-value: 3.698e-05

## 
## Call:
## lm(formula = concern_general ~ (log(MR1 + 10) + log(MR2 + 10) + 
##     log(MR3 + 10) + condition)^2, data = d_moral_merged)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -48.752 -14.595  -4.393  11.230  68.947 
## 
## Coefficients:
##                              Estimate Std. Error t value Pr(>|t|)    
## (Intercept)                  -6093.82    2458.18  -2.479 0.013927 *  
## log(MR1 + 10)                 3021.78    1132.58   2.668 0.008197 ** 
## log(MR2 + 10)                 2713.99     789.50   3.438 0.000702 ***
## log(MR3 + 10)                 -643.70     956.62  -0.673 0.501722    
## conditionrobot                 -26.21      98.37  -0.266 0.790166    
## log(MR1 + 10):log(MR2 + 10)  -1334.25     382.69  -3.486 0.000591 ***
## log(MR1 + 10):log(MR3 + 10)     95.48     286.67   0.333 0.739401    
## log(MR1 + 10):conditionrobot   112.06      69.38   1.615 0.107697    
## log(MR2 + 10):log(MR3 + 10)    183.27     290.03   0.632 0.528104    
## log(MR2 + 10):conditionrobot  -109.43      41.71  -2.624 0.009301 ** 
## log(MR3 + 10):conditionrobot    13.55      36.35   0.373 0.709723    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 21.76 on 220 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.186,  Adjusted R-squared:  0.149 
## F-statistic: 5.028 on 10 and 220 DF,  p-value: 1.384e-06

## 
## Call:
## lm(formula = concern_general ~ (log(MR1 + 10) + log(MR2 + 10) + 
##     log(MR3 + 10))^2 * condition, data = d_moral_merged)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -48.380 -14.461  -4.933  11.436  69.285 
## 
## Coefficients:
##                                            Estimate Std. Error t value
## (Intercept)                                -8278.11    4217.31  -1.963
## log(MR1 + 10)                               3685.53    2034.84   1.811
## log(MR2 + 10)                               3526.34    1242.98   2.837
## log(MR3 + 10)                               -234.07    1525.96  -0.153
## conditionrobot                             -2958.42    4217.31  -0.701
## log(MR1 + 10):log(MR2 + 10)                -1559.42     451.98  -3.450
## log(MR1 + 10):log(MR3 + 10)                   45.01     729.76   0.062
## log(MR2 + 10):log(MR3 + 10)                   52.42     367.99   0.142
## log(MR1 + 10):conditionrobot                 871.50    2034.84   0.428
## log(MR2 + 10):conditionrobot                 935.32    1242.98   0.752
## log(MR3 + 10):conditionrobot                 747.53    1525.96   0.490
## log(MR1 + 10):log(MR2 + 10):conditionrobot  -226.85     451.98  -0.502
## log(MR1 + 10):log(MR3 + 10):conditionrobot   -91.44     729.76  -0.125
## log(MR2 + 10):log(MR3 + 10):conditionrobot  -230.42     367.99  -0.626
##                                            Pr(>|t|)    
## (Intercept)                                0.050937 .  
## log(MR1 + 10)                              0.071491 .  
## log(MR2 + 10)                              0.004985 ** 
## log(MR3 + 10)                              0.878230    
## conditionrobot                             0.483746    
## log(MR1 + 10):log(MR2 + 10)                0.000673 ***
## log(MR1 + 10):log(MR3 + 10)                0.950872    
## log(MR2 + 10):log(MR3 + 10)                0.886867    
## log(MR1 + 10):conditionrobot               0.668865    
## log(MR2 + 10):conditionrobot               0.452579    
## log(MR3 + 10):conditionrobot               0.624719    
## log(MR1 + 10):log(MR2 + 10):conditionrobot 0.616244    
## log(MR1 + 10):log(MR3 + 10):conditionrobot 0.900405    
## log(MR2 + 10):log(MR3 + 10):conditionrobot 0.531875    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 21.86 on 217 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.1894, Adjusted R-squared:  0.1408 
## F-statistic:   3.9 on 13 and 217 DF,  p-value: 1.09e-05

## 
## Call:
## lm(formula = concern_general ~ ((log(MR1 + 10) + log(MR2 + 10) + 
##     log(MR3 + 10))^3) * condition, data = d_moral_merged)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -42.07 -15.15  -4.49  10.45  68.47 
## 
## Coefficients:
##                                                          Estimate
## (Intercept)                                               -202917
## log(MR1 + 10)                                               89596
## log(MR2 + 10)                                               87469
## log(MR3 + 10)                                               84278
## conditionrobot                                            -100257
## log(MR1 + 10):log(MR2 + 10)                                -38599
## log(MR1 + 10):log(MR3 + 10)                                -37247
## log(MR2 + 10):log(MR3 + 10)                                -36390
## log(MR1 + 10):conditionrobot                                45842
## log(MR2 + 10):conditionrobot                                42247
## log(MR3 + 10):conditionrobot                                42505
## log(MR1 + 10):log(MR2 + 10):log(MR3 + 10)                   16076
## log(MR1 + 10):log(MR2 + 10):conditionrobot                 -19331
## log(MR1 + 10):log(MR3 + 10):conditionrobot                 -19411
## log(MR2 + 10):log(MR3 + 10):conditionrobot                 -17951
## log(MR1 + 10):log(MR2 + 10):log(MR3 + 10):conditionrobot     8203
##                                                          Std. Error
## (Intercept)                                                  115903
## log(MR1 + 10)                                                 51666
## log(MR2 + 10)                                                 49811
## log(MR3 + 10)                                                 50196
## conditionrobot                                               115903
## log(MR1 + 10):log(MR2 + 10)                                   22191
## log(MR1 + 10):log(MR3 + 10)                                   22371
## log(MR2 + 10):log(MR3 + 10)                                   21570
## log(MR1 + 10):conditionrobot                                  51666
## log(MR2 + 10):conditionrobot                                  49811
## log(MR3 + 10):conditionrobot                                  50196
## log(MR1 + 10):log(MR2 + 10):log(MR3 + 10)                      9607
## log(MR1 + 10):log(MR2 + 10):conditionrobot                    22191
## log(MR1 + 10):log(MR3 + 10):conditionrobot                    22371
## log(MR2 + 10):log(MR3 + 10):conditionrobot                    21570
## log(MR1 + 10):log(MR2 + 10):log(MR3 + 10):conditionrobot       9607
##                                                          t value Pr(>|t|)
## (Intercept)                                               -1.751   0.0814
## log(MR1 + 10)                                              1.734   0.0843
## log(MR2 + 10)                                              1.756   0.0805
## log(MR3 + 10)                                              1.679   0.0946
## conditionrobot                                            -0.865   0.3880
## log(MR1 + 10):log(MR2 + 10)                               -1.739   0.0834
## log(MR1 + 10):log(MR3 + 10)                               -1.665   0.0974
## log(MR2 + 10):log(MR3 + 10)                               -1.687   0.0930
## log(MR1 + 10):conditionrobot                               0.887   0.3759
## log(MR2 + 10):conditionrobot                               0.848   0.3973
## log(MR3 + 10):conditionrobot                               0.847   0.3981
## log(MR1 + 10):log(MR2 + 10):log(MR3 + 10)                  1.673   0.0957
## log(MR1 + 10):log(MR2 + 10):conditionrobot                -0.871   0.3847
## log(MR1 + 10):log(MR3 + 10):conditionrobot                -0.868   0.3865
## log(MR2 + 10):log(MR3 + 10):conditionrobot                -0.832   0.4062
## log(MR1 + 10):log(MR2 + 10):log(MR3 + 10):conditionrobot   0.854   0.3942
##                                                           
## (Intercept)                                              .
## log(MR1 + 10)                                            .
## log(MR2 + 10)                                            .
## log(MR3 + 10)                                            .
## conditionrobot                                            
## log(MR1 + 10):log(MR2 + 10)                              .
## log(MR1 + 10):log(MR3 + 10)                              .
## log(MR2 + 10):log(MR3 + 10)                              .
## log(MR1 + 10):conditionrobot                              
## log(MR2 + 10):conditionrobot                              
## log(MR3 + 10):conditionrobot                              
## log(MR1 + 10):log(MR2 + 10):log(MR3 + 10)                .
## log(MR1 + 10):log(MR2 + 10):conditionrobot                
## log(MR1 + 10):log(MR3 + 10):conditionrobot                
## log(MR2 + 10):log(MR3 + 10):conditionrobot                
## log(MR1 + 10):log(MR2 + 10):log(MR3 + 10):conditionrobot  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 21.78 on 215 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.2027, Adjusted R-squared:  0.1471 
## F-statistic: 3.644 on 15 and 215 DF,  p-value: 1.145e-05

## 
## Call:
## lm(formula = concern_general ~ ((log(MR1 + 10) + log(MR2 + 10) + 
##     log(MR3 + 10) + condition)^3), data = d_moral_merged)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -40.863 -14.966  -5.071  10.842  68.319 
## 
## Coefficients:
##                                             Estimate Std. Error t value
## (Intercept)                                -122654.7    67756.7  -1.810
## log(MR1 + 10)                                53355.2    29437.6   1.812
## log(MR2 + 10)                                53115.8    29346.4   1.810
## log(MR3 + 10)                                49624.8    29518.7   1.681
## conditionrobot                               -1365.5     4303.7  -0.317
## log(MR1 + 10):log(MR2 + 10)                 -23088.3    12737.1  -1.813
## log(MR1 + 10):log(MR3 + 10)                 -21596.7    12816.5  -1.685
## log(MR1 + 10):conditionrobot                  1764.8     2093.9   0.843
## log(MR2 + 10):log(MR3 + 10)                 -21562.7    12785.4  -1.687
## log(MR2 + 10):conditionrobot                  -265.1     1426.7  -0.186
## log(MR3 + 10):conditionrobot                  -329.9     1647.6  -0.200
## log(MR1 + 10):log(MR2 + 10):log(MR3 + 10)     9379.6     5545.8   1.691
## log(MR1 + 10):log(MR2 + 10):conditionrobot    -387.6      460.0  -0.843
## log(MR1 + 10):log(MR3 + 10):conditionrobot    -320.3      739.1  -0.433
## log(MR2 + 10):log(MR3 + 10):conditionrobot     460.0      548.5   0.839
##                                            Pr(>|t|)  
## (Intercept)                                  0.0717 .
## log(MR1 + 10)                                0.0713 .
## log(MR2 + 10)                                0.0717 .
## log(MR3 + 10)                                0.0942 .
## conditionrobot                               0.7513  
## log(MR1 + 10):log(MR2 + 10)                  0.0713 .
## log(MR1 + 10):log(MR3 + 10)                  0.0934 .
## log(MR1 + 10):conditionrobot                 0.4003  
## log(MR2 + 10):log(MR3 + 10)                  0.0931 .
## log(MR2 + 10):conditionrobot                 0.8528  
## log(MR3 + 10):conditionrobot                 0.8415  
## log(MR1 + 10):log(MR2 + 10):log(MR3 + 10)    0.0922 .
## log(MR1 + 10):log(MR2 + 10):conditionrobot   0.4004  
## log(MR1 + 10):log(MR3 + 10):conditionrobot   0.6652  
## log(MR2 + 10):log(MR3 + 10):conditionrobot   0.4026  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 21.77 on 216 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:    0.2,  Adjusted R-squared:  0.1481 
## F-statistic: 3.857 on 14 and 216 DF,  p-value: 7.403e-06

## Warning: Removed 1 rows containing non-finite values (stat_smooth).

## Warning: Removed 1 rows containing missing values (geom_point).

Start messing with specific moral concern

## Data: d_moral_merged3
## Models:
## r1_spec_raw: concern_score ~ MR1 + MR2 + MR3 + concern_type + condition + 
## r1_spec_raw:     (1 | subid)
## r2_spec_raw: concern_score ~ (MR1 + MR2 + MR3) * concern_type + condition + 
## r2_spec_raw:     (1 | subid)
## r3_spec_raw: concern_score ~ (MR1 + MR2 + MR3 + concern_type)^2 + condition + 
## r3_spec_raw:     (1 | subid)
## r4_spec_raw: concern_score ~ (MR1 + MR2 + MR3 + concern_type + condition)^2 + 
## r4_spec_raw:     (1 | subid)
## r5_spec_raw: concern_score ~ (MR1 + MR2 + MR3 + concern_type)^2 * condition + 
## r5_spec_raw:     (1 | subid)
## r6_spec_raw: concern_score ~ (MR1 + MR2 + MR3 + concern_type + condition)^3 + 
## r6_spec_raw:     (1 | subid)
## r7_spec_raw: concern_score ~ (MR1 + MR2 + MR3 + concern_type)^3 * condition + 
## r7_spec_raw:     (1 | subid)
##             Df   AIC   BIC logLik deviance    Chisq Chi Df Pr(>Chisq)    
## r1_spec_raw 10 22291 22349 -11136    22271                               
## r2_spec_raw 19 21998 22109 -10980    21960 310.7760      9  < 2.2e-16 ***
## r3_spec_raw 22 22003 22131 -10979    21959   1.4891      3  0.6847954    
## r4_spec_raw 28 22002 22165 -10973    21946  12.7688      6  0.0468574 *  
## r5_spec_raw 40 22000 22233 -10960    21920  26.1706     12  0.0101527 *  
## r6_spec_raw 50 21986 22278 -10943    21886  33.4061     10  0.0002327 ***
## r7_spec_raw 60 21993 22342 -10936    21873  13.7322     10  0.1855586    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## Linear mixed model fit by REML ['lmerMod']
## Formula: concern_score ~ (MR1 + MR2 + MR3 + concern_type + condition)^3 +  
##     (1 | subid)
##    Data: d_moral_merged3
## 
## REML criterion at convergence: 21782.1
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -3.2972 -0.5626  0.0051  0.4418  3.9876 
## 
## Random effects:
##  Groups   Name        Variance Std.Dev.
##  subid    (Intercept) 323.3    17.98   
##  Residual             283.0    16.82   
## Number of obs: 2515, groups:  subid, 231
## 
## Fixed effects:
##                                    Estimate Std. Error t value
## (Intercept)                         26.4008     6.7751   3.897
## MR1                                 22.8547     7.1126   3.213
## MR2                                  2.9675     3.4237   0.867
## MR3                                 -8.8430     6.3475  -1.393
## concern_typephy                      8.0218     2.6136   3.069
## concern_typesoc                     -2.4300     2.6127  -0.930
## concern_typeper                     -4.6537     2.6226  -1.774
## conditionrobot                      18.2880     6.7818   2.697
## MR1:MR2                             -4.4060     3.4802  -1.266
## MR1:MR3                              0.3335     6.8014   0.049
## MR1:concern_typephy                  1.2414     2.6882   0.462
## MR1:concern_typesoc                 -5.9015     2.6881  -2.195
## MR1:concern_typeper                 -2.2884     2.6973  -0.848
## MR1:conditionrobot                   4.9946     7.1086   0.703
## MR2:MR3                              0.3147     3.0035   0.105
## MR2:concern_typephy                  0.3706     0.9983   0.371
## MR2:concern_typesoc                  1.0892     0.9979   1.091
## MR2:concern_typeper                 -0.9955     1.0014  -0.994
## MR2:conditionrobot                  -2.4499     3.3787  -0.725
## MR3:concern_typephy                  1.5569     1.1110   1.401
## MR3:concern_typesoc                  0.4154     1.1101   0.374
## MR3:concern_typeper                  0.4242     1.1152   0.380
## MR3:conditionrobot                   3.2067     6.3625   0.504
## concern_typephy:conditionrobot      -5.4281     2.6022  -2.086
## concern_typesoc:conditionrobot       0.3782     2.6033   0.145
## concern_typeper:conditionrobot      -3.5330     2.6172  -1.350
## MR1:MR2:MR3                          8.2230     4.3190   1.904
## MR1:MR2:concern_typephy              2.2813     1.3234   1.724
## MR1:MR2:concern_typesoc             -2.7981     1.3235  -2.114
## MR1:MR2:concern_typeper              3.2123     1.3249   2.425
## MR1:MR2:conditionrobot              -5.3550     3.5250  -1.519
## MR1:MR3:concern_typephy              1.4061     1.1705   1.201
## MR1:MR3:concern_typesoc             -3.0578     1.1703  -2.613
## MR1:MR3:concern_typeper              1.6078     1.1740   1.369
## MR1:MR3:conditionrobot              -6.4390     6.8990  -0.933
## MR1:concern_typephy:conditionrobot  -8.9588     2.7648  -3.240
## MR1:concern_typesoc:conditionrobot   5.9755     2.7629   2.163
## MR1:concern_typeper:conditionrobot  -1.6125     2.7692  -0.582
## MR2:MR3:concern_typephy              1.2545     1.0674   1.175
## MR2:MR3:concern_typesoc              0.5214     1.0672   0.489
## MR2:MR3:concern_typeper             -0.1657     1.0679  -0.155
## MR2:MR3:conditionrobot               8.3439     4.3498   1.918
## MR2:concern_typephy:conditionrobot   5.0086     1.4704   3.406
## MR2:concern_typesoc:conditionrobot  -2.0552     1.4700  -1.398
## MR2:concern_typeper:conditionrobot   0.9625     1.4734   0.653
## MR3:concern_typephy:conditionrobot   0.3889     1.5321   0.254
## MR3:concern_typesoc:conditionrobot  -2.1523     1.5307  -1.406
## MR3:concern_typeper:conditionrobot  -0.5007     1.5426  -0.325

Plots with specific moral concern

## Warning: Removed 37 rows containing non-finite values (stat_smooth).

## Warning: Removed 37 rows containing missing values (geom_point).

## Warning: Removed 37 rows containing non-finite values (stat_smooth).

## Warning: Removed 37 rows containing missing values (geom_point).

## Warning: Removed 37 rows containing non-finite values (stat_smooth).

## Warning: Removed 37 rows containing missing values (geom_point).

## Warning: Removed 26 rows containing non-finite values (stat_smooth).

## Warning: Removed 26 rows containing missing values (geom_point).

## Warning: Removed 26 rows containing non-finite values (stat_smooth).

## Warning: Removed 26 rows containing missing values (geom_point).

## Warning: Removed 26 rows containing non-finite values (stat_smooth).

## Warning: Removed 26 rows containing missing values (geom_point).